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Q^- 1 A new procedure for solving linear differential equations 

*0 ■ Linear differential equations play a crucial role in various branches of science and mathematics. 

Second order differential equations routinely manifest in the study of quantum mechanics, in 
\ connection with Schrodinger equation. There are various techniques available to solve a given 

' differential equation, e.g., power series method, Laplace transforms, etc. Not many general 

■ methods applicable to differential equations of arbitrary order exist in the literature (see [1] and 

references therein). We make use of a general method for solving linear differential equations 
of arbitrary order to construct new representations for the solutions of the known second order 
linear differential equations, both without and with a source term. This method has found 
applications in solving linear single and multi variable differential equations. The solutions 
of linear differential equation with a source and development of a new Quasi-Exactly Soluble 
(QES) system are the new results of this paper. 

1.1 Case (i). Linear differential equation without a source term 

After appropriate manipulation, any single variable linear differential equation can be brought 
to the following form 

[F{D)+P{x,d/dx)]y{x)=Q, (1) 

where D = x-^, = Xln '^"-^'^ ^ diagonal operator in the space of monomials spanned 

by and OnS are some parameters. Here P{x, d/dx) = ^ ■ • Ci,jX^{-^y , where Cij = if i = j. 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue 'Integrable Systems and Related Topics'. 
The full collection is available at |http: / /www. emis.de/journals/SIGM A /kuznetsov. html. 
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Notice that, since F{D) is a diagonal operator, 1/F{D) is also well defined in the space of 
monomials. The following ansatz 



1 



F{D) 



P{x, d/dx) 



= CxGxx^ 



(2) 



' dx' 



X 



is a solution of the above equation, provided F{D)x^ = and the coefficient of x'^ in P{x 
should be zero [6]. In order to realise ([2]) as a power series we impose the requirement that 
P{x, ^) lowers (i.e. Cj.j = for i > j) or raises the degree of monomials. Substituting equa- 
tion ([2]), modulo Ca, in equation ([T]) 
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x^ = 0. 



Equation ^ connects the solution of a given differential equation to the monomials. In 
order to show that, this rather straightforward procedure indeed yields non-trivial results, we 
explicitly work out a few examples. Consider the Hermite differential equation, which arises in 
the context of quantum harmonic oscillator. 



1 

^-^-2d^ 



Hn{x) = 0. 



Here F{D) = D — n and F{D)x'^ = yields X = n. Hence 



m=0 



1 



D-n 



{-l/2){d^/dx'^) 



X . 



Using [D, {d^ /dx^)] = —2{d^/dx^) it is easy to see that 



1 



[D-n] 



-l/2){d^/dx^) 



m -. 

{-l/2)"^{d'/dx')^l[—-x- 
1=1 ^ 



and 



oo 

Hn{x) = C„ V(-l/4r — (dVdx2)™x" = One 
' m! 



1 d-^ 

4 d^x", 



m=0 



this is a well-known result. Similar expression also holds for the Lagurre polynomial which 
matches with the one found in [5]. Below, we list new representations for the solutions of some 
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frequently encountered differential equations in various branches of physics and mathematics [1] . 
Notice that, the cases of confluent hyper geometric, hypergeometric, Chebyshev type II and 
Jacobi solutions are given in [6], and reproduced here for the sake of completeness. 
Legendre polynomial 

p^(^) ^ (^^^g-{l/(2[D+n+l])}{d2/dx2)^n_ 

Associated Legendre polynomial 

P^i^x) = C„(l - 3;2)nV2g-{l/(2[Z)+n+m+l])}(d2/'^a:2)^n-m_ 

Bessel function 

Generalized Bessel function 

Gegenhauer polynomial 

Hypergeometric function 

y±(a,/?;7;x)=C±e-«V(^+^±)>V^^ 

where \± is either a or (3 and A = x-^ + 7^- All the above series solutions have descending 
powers of x. In order to get the series in the ascending powers, one has to replace x by ^ in the 
original differential equation and generate the solutions via equation ([2]). However, the number 
of solutions will remain the same. One can also generate the series solutions by multiplying the 
original differential equations with x^, and then, rewriting x^-^ = D{D — 1) = F{D). 
The solution for the following equation with periodic potential 

^-^ + acos(x)y = (3) 

can be found after multiplying equation ([3]) by and rewriting x"^-^ as {D — 1)D to be 



E ^ n 



^-^ ml I (2nj)\ 

i,{n,}=0 U=l ^ 

m 

n 







m+l—r 

m + A/2 - r + Yl '^i 

i=l 






m+l—r 

m + A/2 + 1 - r + 

i=l 


} 



, m . 

X '=1 



where A = or 1. In the same manner, one can write down the solutions for the Mathieu's 
equation as well. 

Chebyshev polynomials Tn{x) and Un{x) 



r 1 1 dM 

Tn{x) = Cue I2(l3+n)d.2/a;n 
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and 



f 1 1 d2 



where Cn and are appropriate normalization constants. 
Jacobi polynomial 



m=0 



1 d_ 

{D-n){D + a + f3 + l)\dx'^^^^ "^'dx 



Schlafli, Whittaker and for that matter, any solution of a second order linear differential equation 
without a source term can be solved in a manner identical to the above cases. 



1.2 Case (ii). Linear differential equation with a source term 

Consider, 

{F{D) + P{x, d/dx)) y{x) = Q{x). 
Now the solution can be found in a straightforward way as 
1 



P{x,d/dx) y{x) = Q{x), 



and 



OO 

E 

m=0 



i-iy 



[F{D) 



P{x, d/dx) 



F{D) 



Q{x). 



We list below a few examples for the above case. 
Neumann's polynomial 



OJx) 



.r=0 



1 



[(D + l) 



tX 



[{D + 1)2 - n2] 



X (xcos^(n7r/2) + nsin^(n7r/2)). 
Lommel function 

OO 

w{x) = Y,{-l) 



m=0 



(D2 _ „2 



The cases of Struve, Anger and Weber functions are identical to the above ones. 

It is a priori not transparent that all the solutions of a given differential equation can be 
obtained through the present approach. If F is a polynomial of the same degree as the order 
of the differential equation, and has distinct roots, the linearly independent solutions can be 
obtained through this approach. As has been pointed out in an earlier paper [8] a given single 
variable differential equation can be cast in the desired form, as required by the present approach, 
in more than one way through multiplication by powers of x. In a number of cases, these lead 
to different solutions. However, the case of degenerate roots, as also the case of inhomogeneous 
equations needs separate consideration, which we hope to investigate in future. 
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2 A new quasi-exactly solvable model 

The above procedure is also applicable to differential equations having higher number of sin- 
gularities e.g., Heun equation and its generalization [7]. The same can be used to generate 
quasi-solvable models. It is worth mentioning a quasi-solvable model based on Heun's equation 
has been studied; it has been shown that, this system lacks the SL{2, R) of many well-known 
QES model [U [9] . We now proceed to generalize this system for obtaining a new QES potential. 
The differential equation under consideration is given by: 

r + i±f + 'I /- . -% f = 0. (4) 

{2x x — l 2[x + e-^) ] x[x — l){x + e^) 
When = the above reduces to Heun equation. Under the following change of variable 
sinh^ f 



1 + 1 + sinh2 ' 



and with ^' = (1 — xY f[x) the above equation can be cast in the form of Schrodinger eigenvalue 
problem. Here s = (1 — E/ p^Yl"^^ where E is energy of the system. The constants a, /? and q 
in the equation ^ are related to s in the following form: 

a = -^--s, P=\-^ and q = {I - s^){l + e^) - ]^se^ - -^{l - 2e^) . 
The corresponding potential is given by 

J 8sinh^ f - - 1) sinh^ f + 2(^ - ^ - 2) ] 



8(1 + ^ + sinh2 f )2 e2sinh2f j" 

For the purpose of finding solutions we cast the generalized Heun equation Q in the form, 
[F{D) + P]y{x) = 0. Multiplying by x^, we get 



_4xV-^-2e2x-^-40 
dx^ dx 



4x3(e2 - 1)-^ + 2(3e2 - 2 + 4se^)x''^ 
dx'' dx 



4qx + 4a;'^— r + 8(1 + s)x^— + 4a/3x^ 
dx"' dx 



fix) = 0. 

In the above equation F{D) = —Ae^D"^ + 2e^D — Aft and the condition F{D)x^ = yields 



e± = ^ (i ± ^1 - wn/e^ 

For polynomial solutions either of the roots of the equation must be an integer. Taking ^_ = m 
we obtain the following relation for the allowed values of m 

2m2 - m + 2n/e^ = 0. 

Let us consider a case in which n = -eV2. Taking care of normalization condition and 
positivity of energy, the allowed value of m is 1. This yields a solution of the form: 

f{x) = aix + a2X^ -|- a^x^ + • • • . 

If the polynomial terminates at x^~^ then the coefficients a„ and a„+i should be zero. Using 
this constraint, we get the following condition for n 

+ i2n -l) + n^ -n-^ = 0, 



6 



N. Gurappa, P.K. Jha and P.K. Panigrahi 



giving the allowed values of n as 2, 3. For n = 2, s = l/2, i? = and = 1 we get the wave 



4 

function 



^(y) = N 



2 / smh'{py/2) \ 



2 + sinh2(py/2) / V2 + smh'^{py/2) J ' 



where N is the normalization constant. For n = 3, s = — ^, which makes the wave function 
un-normalizable. This procedure may find applicability to unravel the symmetry properties of 
Heun equation and its generalization [2]. We have earlier studied the symmetry properties of 
confluent Hypergeometric and Hypergeometric equations using the present approach [8]. The 
fact that, the solutions are connected with monomials makes the search for the symmetry rather 
straightforward. 

In the case when 7^ 0, in equation the additional term in Heun equation manifests in 
the potential in the Schrodinger equation, as a term having ^ type singularity at the origin. 
This is very interesting since it is of the Calogero-Sutherland type, a system much investigated 
in the literature. Like the Calogero-Sutherland case, a Jastrow type factor in the wave function 
appears because of this singular interaction. At a formal level, this interaction modifies the 
singularity at the origin without adding any new singularity. We note that x = is still 
a regular singular point, since the limit x^-^ as x ^ is finite. 

In conclusion, we have developed a new method to solve linear differential equations of 
arbitrary order, both without and with source term. Using this method, we worked out a few 
examples for the case of second order linear differential equations. We obtained known, as 
well as, the new representations for the corresponding solutions. The same approach is also 
applicable to QES models based on Heun equation and its generalizations. In particular a new 
potential of QES type is constructed through this approach. We intend to analyze the symmetry 
properties of the Heun equation through the present formalism. 
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